Energy-gap analysis of quantum spin-glass 
transitions at zero temperature 
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Abstract. We study for random quantum spin systems the energy gap between the ground 
and first excited states to clarify a relation to the spin-glass-paramagnetic phase transition. We 
£h , find that for the transverse Sherrington-Kirkpatrick model the vanishing of the averaged gap 

Ch ■ does not identify the transition. A power-law form of the gap distribution function leads to 

iy>j ' power-law distributions of the linear, spin-glass, and nonlinear susceptibilities. They diverge 

at different points, which we attribute to a quantum "Griffiths" mechanism. On the other 
hand, the classical treatment is justified for the transverse random energy model and the phase 
transition can be found by a sudden change of the ground state. 

a' 

a . 

1. Introduction 

Spin glasses are known as nontrivial states of random spin systems and a lot of studies have 
revealed their interesting properties [HO [3]. Now, the mean field theory of the classical spin-glass 
has been established and is recognized as a useful tool for attacking various random systems. 
As one of such applications, quantum systems have been studied for many years since the work 
C\| \ by Bray and Moore [I]. Recently, much attention has been paid to quantum systems for their 

applications to information processing such as quantum annealing [5j [6j [7]. In contrast to 
classical systems, mean field analysis has a serious limitation and a reliable analytical method 



is needed. 



At finite temperature, state of the system is determined by the free energy. A competition 
between energy and entropy induces a phase transition between different phases. At zero 
temperature, the system is determined by the energy only and the entropy does not play any 
role. The quantum fluctuations become important at low temperature and the system falls 
into its ground state. Then, the phase transition is realized by a change of the ground state, 
which can be found by looking at the energy gap between the ground and first excited states. 
Although the mechanism of a quantum phase transition is simpler than that of a thermodynamic 
phase transition, it is rather difficult to treat quantum fluctuation effects even for clean systems 
[8]. The effect of disorder further makes the analysis even more difficult. As in the classical 
system, the mean field theory of quantum spin-glass systems has been developed. However, a 
lot of studies use the static approximation [1] which neglects the quantum fluctuation effects. 
Although this approximation is justified for not-too-low temperatures, no reliable result can be 
obtained at zero temperature. 

In the present paper, we give a detailed study of quantum phase transitions in terms of the 
energy gap |2J. As typical random quantum systems which exhibit a spin-glass phase transition, 
we treat the Sherrington-Kirkpatrick (SK) model [TU] and random energy model (REM) [11 1 



with a transverse field. The transverse REM can be solved exactly and shows a discontinuous 
transition [12] . The transverse SK model has not been solved exactly, but many studies indicate 
a continuous transition. Since these models without transverse field can be treated exactly by 
a mean field theory, most of previous studies took a similar approach. We study these systems 
by using the energy gap analysis. In section [2j we define the models and consider the energy 
gap by using perturbation theory. The energy gap is also examined numerically and we see 
that the gap strongly fluctuates from sample to sample. In order to see a phase transition, 
several kinds of susceptibilities are studied in terms of the energy gap in section [3J We find 
there a strong quantum fluctuation effect. After discussing a sufficient condition for this effect 
to occur in section |U we discuss its mechanism in section [SJ The final section O is devoted to 
give conclusions, perspectives, and possible applications. 

2. Energy level distribution 

In this section, we define quantum spin-glass models and examine the energy gap between the 
ground and first excited states. 

2.1. Model 

We treat the p-body interacting spin-glass model with transverse field 

N N 

H = ~ £ Jhi^i p a! 1 al---af p -TY / of, (1) 

*i<«2<---<*p i=l 

where a\ ,x are Pauli matrices on site i, N is the number of the site, and T is the transverse field. 
We take an ensemble average over random interaction Ji 1 i 2 ...j J , with a probability distribution 



Np- 1 I Nf- 1 , 2 | 
P(^i-v) = y^yj2exp|--^(J ili2 ... ip ) j, (2) 

where J is the strength of the random fluctuations. Here we study the cases of p = 2 (SK model) 
and p — > oo (REM), and the case of other-p will be reported elsewhere. The quantum REM can 
be solved exactly and at zero temperature the system shows a discontinuous phase transition 
between the spin-glass and paramagnetic phases at T/J = Vhi2 ~ 0.83 |12j . The SK model has 
not been solved exactly, but a lot of studies reported a continuous transition around the point 
T/J = 1 ~ 2 when the temperature is equal to zero. A naive perturbative calculation gives 
T/J = 1 and the static approximation gives T/J = 2 |13|. 114] . while more advanced methods 
indicate a point around T/J ~ 1.5 [EJ El OH [TBI Q21 [20]. 

2. 2. Perturbative calculation of the energy gap 

Here we estimate for the SK model the energy gap between the ground and first excited states 
by using a perturbative expansion with respect to J/T. When the randomness J is absent, the 
Hamiltonian is exactly solved to give energy levels E^ = —(N — 2n)T with nonnegative integer 
n. Each of the levels has a degeneracy N\/n\(N — n)\ which is split into many nondegenerate 
ones for finite- J. The unperturbed ground state is denoted by n = and the first excited 
state belongs to the sector with n = 1. The Hamiltonian of the n = 1 sector is given by 
HQ = —(N — 2)r — Jij. Since Jy are Gaussian random variables, the energy level distribution 
is given by that of the Gaussian orthogonal ensemble [21] . The energy levels distribute in a 
semicircle form and the lowest energy level is given by E\ = —(N — 2)T — 2 J. Neglecting the 
correlations between blocks with different n's, we obtain the energy gap A ~ 2(r — J) and the 
degenerate point T = J. 



The subleading corrections are calculated as follows. The ground state n = is coupled 
to N(N — l)/2-levels in the sector n = 2 by random interactions Jij. Up to second order in 
perturbation theory, we obtain in average, 

In the same way, the first excited state is directly coupled to the sector n = 3. The correction 
is roughly estimated as 

E f) „ _ {N _ 2 )r - 2J - £(N - 3). (4) 
As a result we have the energy gap 

A~2(r-J) + ^, (5) 

and the degenerate point T/J = (2 + \/2)/4 ~ 0.85. Although the present analysis is based on a 
perturbative expansion and the degenerate point cannot be estimated exactly, it is probable to 
say that the averaged degenerate point does not coincide with the phase transition point around 
1.5. On the other hand, a similar perturbative calculation was carried out for the transverse 
REM [22J and a good estimate of the phase transition point was obtained. 



2. 3. Numerical result of the energy gap 

We estimate the energy gap by using the numerical calculation. The Hamiltonian of the size 
2^ is diagonalized by using the Lanczos method and the average is taken over more than 20000 
samples. The number of the site is taken to be N = 8 — 22. We note that for the SK model, 
the Hamiltonian commutes with P = J\i=i a f an d the matrix splits into two blocks. This allows 
us to diagonalize their two matrices separately. For the REM, we take the Derrida's original 
Hamiltonian for the random part and do not use the original Hamiltonian ([1]) with p — > oo. 
The Gaussian random variables with the average [Ej\ = and the variance [Ef] = NJ 2 /2 are 
produced 2^ times and are used for the diagonal part of the Hamiltonian in ^-basis. 

In figure [H we plot the averaged energy gap. For the SK model, the gap is an increasing 
function of T and becomes negligible at values smaller than T/J ~ 0.6. This result is consistent 
with the perturbative calculation. From this figure, we see that it is hard to identify a phase 
transition point around T/J ~ 1.5 even if the fluctuations are taken into account. On the other 
hand, the gap shows a complicated behavior for the REM. The averaged gap becomes minimum 
around the phase transition point and the fluctuations behave differently on both sides of the 
point. These results imply that it is possible to find a phase transition from the averaged energy 
gap for the REM, while it is hard for the SK model. 

In order to see the difference more in detail, we show the energy gap distribution function 
in figure [2j For the SK model, we find long tails of the distribution function. It is hard to 
imagine that these fluctuations vanish at the thermodynamic limit and we expect a significant 
role of fluctuations. A more complicated behavior can be observed for the REM. A single peak is 
replaced by another one, which implies a discontinuous transition between two different states. 

These analytical and numerical results show that for the SK model the fluctuations play a 
significant role and the averaged gap is not a useful quantity to describe the phase transition. 
Concerning the REM, the averaged gap is expected to be useful although the fluctuations are 
large as well. Therefore, it is important to understand where this difference comes from. 
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Figure 1. Averaged energy gap A = [E\ — 
Eq\ between the ground and first excited 
states for the SK model and REM with 
N = 22. The variance of the energy gap 
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Figure 2. Gap distribution function for the 
SK model and REM with N = 22. For the 
REM, T/J ~ 0.83 corresponds to the phase 
transition point. 



3. Susceptibility 

In the standard theory of critical phenomena, phase transitions are identified by a singularity 
of response functions such as the magnetic susceptibility. In this section, we discuss the linear, 
spin-glass, and nonlinear susceptibilities. For our purpose to find a transition in terms of the 
energy gap, the spectral representation is useful and is discussed in detail. 



3.1. Spectral representation 

The magnetization of the system is given by the response of the longitudinal magnetic field h 
whose contribution to the Hamiltonian is given as H —> H — h J2i °f • We have 

where Z = Tre - ^ is the partition function and /3 = 1/T is the inverse temperature. The 
magnetization can be expanded with respect to h as m = x\h — Xn\h 3 + ■ • "j which defines the 
linear (xi) and nonlinear (xni) susceptibilities [231, [23]. 

Our analysis is based on the imaginary time formalism [51 [2S] . The linear susceptibility is 
expressed as 

N 



Xl 



— V 



dr / dr' 
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- Tr e^ H a?(r)a](r') 



Pq, 



(7) 



where af(r) 



HT afe~ HT , 



t is the imaginary time running from to /3, and q is the spin- 



glass order parameter. Inserting the completeness relation of eigenstates of the Hamiltonian, we 
obtain that the first term of the right hand side of equation ([7]) is equal to 
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(8) 



where \n) is an eigenstate with the energy E n . The first term comes from degenerate states and 
the second from nondegenerate ones. The first term is equal to the static part oj = of the 
spectral function 
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For the static part x = x( u = 0), the energy is conserved through a virtual process in imaginary 
time and a classical treatment is possible. Then, the linear susceptibility is expressed by a sum 
of classical and quantum parts. The classical part Xi is equal to — q) and can be treated 
by the static approximation. It corresponds to the linear susceptibility for the classical system 
given by x\ = — q) [2]< On the other hand, the quantum part is expressed as the second 
term of equation ([5]). At the zero temperature limit, we assume that the ground state |0) with 
the energy Eq is nondegenerate and obtain for the quantum part 



(q) 
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E n — En 



(10) 



In the same way, the quantum part of the nonlinear susceptibility is given by 
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(11) 



In classical spin-glass systems, the nonlinear susceptibility is an important quantity since it 
is related directly to the spin-glass susceptibility Xsg as Xnl = /3(Xsg — 2/3 2 /3). A spin-glass 
transition which is characterized by the divergence of Xsg can be observed directly by %n\- 
However, in the present case, the quantum part of the spin-glass susceptibility is given by 



y (q) 
Asg 
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l(0|E,<|n)|' 
E n — En 



(12) 



X^x and xig^ behave differently and there is no direct relation between them. The spin-glass 
susceptibility is not equivalent to the nonlinear susceptibility, which is a crucial difference from 
the classical system. 

If the transition is of quantum nature, the singularity comes from the quantum part. We see 
from above expressions that the three susceptibilities are divergent when the gap between the 
ground and first excited states collapses into zero. This is understood as the standard mechanism 
of a quantum phase transition as we mention in section [TJ However, in random systems, the gap 
vanishing point fluctuates from sample to sample. If the fluctuation is small, the transition point 
can be given by the averaged gap vanishing point. But, as we see in section [21 the fluctuation 
is not negligible and the vanishing point does not coincide with the transition point reported in 
previous studies. If the gap distribution function behaves in a power-law form P(A) ~ A fc , we 
can find that three susceptibilities diverge at different points. Since xu Xsg, and Xni are roughly 
proportional to the gap as 1/A 1,2 ' 3 respectively, they diverge at k = 0, 1, 2 respectively. We note 
that this behavior is specific to random and quantum systems. The quantum nature is necessary 
to obtain the spectral representations and a power-law distribution of the gap can be obtained 
only for random systems. 
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Figure 3. Inverse susceptibilities (quantum 
part) for the SK model with N = 18. 
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Figure 4. The cumulative distribution func- 
tion of the energy gap F(A) = J A dA'P(A') 
(N = 12,14,16,18,20, and 22 from below to 
above) for the SK model. 



3.2. Numerical result 

In order to see whether the above scenario holds, we analyze numerically the spectral 
representation of the susceptibilities. In figure El we plot the quantum part of \h Xsg-> an d 
Xni f° r the SK model. They become very large at small transverse fields and there is a relevant 
difference between Xsg an d X$ ■> which support our expectation. We also plot the cumulative 
distribution function of the energy gap in figure [H The size dependence of the function implies 
that the power-law persists in the thermodynamic limit. 

We also plot the distribution function of the susceptibilities in figure [5j The distribution 
function P(x) shows a power-law dependence at large-x as P(x) ~ X~ B • 111 this form of 
the distribution function, the averaged susceptibility diverges when s < 1. We analyze the 
distribution functions for each of the susceptibilities and find that the divergent point can be 
different with each other, which is consistent with the analysis of the gap distribution. Although 
our system size is not so large enough to give a reliable extrapolation value, our result of the spin- 
glass transition, identified by the divergence of Xsg> supports the results in pT5l [TBI fTTl fTHl fT9l. [20] . 
Previous studies did not pay any attention to the difference between the spin-glass and nonlinear 
susceptibilities and it may be worth reconsidering the results more carefully. 

We note that the present numerical analysis is only for the quantum part. To find an 
observable physical value of susceptibilities, we must include the classical part. The classical part 
of the susceptibilities can be analyzed by the static approximation. For the linear susceptibility, 
= /3(x—q) on ly gives a cusp and does not show any strong singularity. The classical nonlinear 
and spin-glass susceptibilities can be divergent and we must study closely which part between 
the classical and quantum parts is important. It can be considered that the classical part is 
important at the spin-glass phase and the quantum part is dominant at the paramagnetic phase. 
We treat a related problem in the next section and a detailed study is presented elsewhere. 

3.3. Effect of higher excited states 

Up to here, our attention has been paid to the gap between the ground and first excited states. 
The spin-glass system has a complicated ground state which is described by a replica-symmetry- 
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Figure 5. Susceptibility distribution func- 
tions for the SK model with N = 18. The 
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Figure 6. Spectral function x( w ) f° r the 
SK model (N = 20). The shaded parts are 
contributions from the first excited state. 



broken order parameter, and it might be natural to expect that at the transition point a vast 
amount of degeneracy occurs at the same time. In order to see effects of higher order excited 
states, we plot the spectral function ([9]) for the SK model in figure El For large-r, the function 
was analytically calculated as [26] 

In this case, the contribution from the first excited state cannot be distinguished from those of 
other states. This property drastically changes when T is small. We see that a single peak from 
the first excited states is separated from other contributions and plays a dominant role for the 
phase transition. Actually, when the transverse field is absent, we have a symmetry that the 
Hamiltonian is unchanged under the gauge transformation of —¥ —of. As a result, the energy 
level is doubly degenerate in the classical limit. Thus, for the present model, the first excited 
state plays a dominant role for the phase transition. 

This result implies that the second term of the right hand side of equation (jlip gives a 
subleading correction since the expression includes matrix elements of higher order excited states 
(|n')). Numerically we confirm that this term is not important to find a divergent point, which 
makes numerical calculations easier. 



(13) 



4. Localization of state 

In the previous section, we have treated the SK model only. The spectral representation is a 
general relation and is independent of the explicit form of the Hamiltonian. However, it is evident 
from the known exact solution that the above scenario cannot be applied to the REM. Actually 
the classical treatment is justified for the REM and the quantum part of the susceptibility 
does not give any significant contributions. This can be understood by examining the inverse 



participation ratio (IPR) denned as 



£lH^|o>r 



(14) 



This function describe how many states are included in the state of|0) = Xm=i c «l n )- ^ ^ 
includes M-states, c n ~ 1/\[M and we obtain / ~ 1/M. Thus this function is a useful quantity 
to see whether the state is localized in eigenstate space |27| . 
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Figure 7. Inverse participation ratio 
(N=12). The summation over n in equation 
()14p is restricted to n = 0(1) for a plot with 
caption n = 0(1). 
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inverse participation ratio (N=12). 



In figure [3 we plot the averaged IPR for both models and find a significant difference between 
them. For the SK model, the IPR decreases uniformly from 1 to 1/N as a function of T. This 
behavior reflects the fact that the gauge symmetry is important at small-r and, at large-r, the 
one spin flip changes the ground state to the sector n = 1 with iV-levels. af |0) mainly consists of 
the first excited state. There is no such symmetry for the REM and the IPR becomes minimum 
around the phase transition point. At large-r the IPR gradually approaches to unity. of|0) 
almost stays in the ground state |0) and the excitations to the higher energy levels are small. 
This result shows that the quantum fluctuation effect can be negligible for the REM. This is the 
reason why the discussion in the previous section cannot be applied to the REM. The state of 
the system is determined by the ground state. 

We also find for the IPR distribution function in figure [8] that the IPR is always larger than 
1/N. This shows that all states 2 N do not contribute to determine the state of the system, 
which may allow us to construct a reduced effective model. As in the energy gap distribution, 
a double peak structure is observed for the REM, while it is not for the SK model. 



5. Quantum Griffiths singularity 

The power-law distribution of the susceptibility was also found in finite-dimensional systems 
[23 ES EQj EU [32j [33]. In two and three dimensional systems, the nonlinear susceptibility 
diverges at a point of T which is larger than the phase transition point. This was understood 



the Griffiths singularity effect [34|. [35] . Rare configurations of ordered clusters in a disordered 
state give a singular contribution to the susceptibility. A similar effect is known in disordered 
electron systems as anomalously localized states [36J. 

Since the origin of the singularity comes from ordered clusters in spatial dimensions, this 
effect should be absent in the present infinite-ranged models. The authors in |29[ 130} l32j [33] 
called this mechanism as the quantum Griffiths singularity. However, the essential mechanism is 
the same as the classical one. The quantum nature is relevant for how this mechanism is related 
to physical quantities. 

Our finding is that, even for infinite-ranged systems, a similar power-law behavior can be 
observed for the susceptibility distribution. We cannot attribute this behavior to the ordered 
clusters as in the finite dimensional case. The effect must be a purely quantum mechanical one 
which can be expressed by fluctuations in imaginary time. Therefore, it is natural to expect that 
ordered clusters are formed in the imaginary time direction. Compared to the spatial dimension, 
the crucial difference is that there is no disorder in this direction. Theoretically, the imaginary 
time correlation arises only when the ensemble average is carried out for the free energy, which 
implies that the present Griffiths-like effect can be obtained for random quantum systems. The 
randomness induces quantum fluctuations, which allows a formation of ordered clusters. If 
we accept this scenario, it is not difficult to obtain a power-law form of the gap distribution 
function. Referred to the classical Griffiths mechanism [8], we expect an exponentially small 
gap A ~ Aoe^ CT /( fc+1 ) in r is formed with a very small probability e~ CT , where Ao, k, and c are 
some constants. Summing over r, we obtain 



Since < r < /3, the power law can be obtained only for large-/?. This feature is consistent with 
our expectation that the present effect is a purely quantum one and becomes important at low 
temperatures. 

This mechanism is not specific to the present infinite-ranged models and the finite dimensional 
systems can share the same property. For the finite dimensional systems, the classical and 
quantum Griffiths singularities appear at the same time. They must be treated in a different 
way and their contributions to physical quantities are expected to be different. It is an important 
issue to be clarified in future studies. 

6. Conclusions 

In conclusion, we have discussed quantum spin-glass phase transitions in terms of the energy 
gap. We found for the transverse SK model that the gap distribution shows a power-law behavior 
near the origin. As a result, depending on its exponent, the linear, spin-glass, and nonlinear 
susceptibilities diverge at different points. On the other hand, the transverse REM does not 
show such behavior. Although the gap distribution shows a broad form, the state of the system 
can be described by the ground state and the classical treatment is justified for the REM. The 
difference can be understood by the symmetry of the system. The role of symmetry can be more 
closely studied by generalizing the model ([T|) to arbitrary-p. Details are presented elsewhere. 

The present results are considered a universal feature that can be seen in broad random 
quantum systems since the mechanism of the quantum fluctuation effects induced by disorder 
is independent of the details of the Hamiltonian. Our result for the SK model shows that the 
self- averaging property does not hold at zero temperature and the quantum phase transition is 
determined by rare configurations of disorder. Random quantum systems can be totally different 
from classical thermodynamic systems. In order to reveal such properties, we consider that the 
present approach is useful and can be a supplement to the dynamical mean field theory as 
discussed in [T71 [20] . Work in this direction is under progress and we hope that such universal 
properties are clarified in near future. 
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